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Abstract 

A 1-meter-long trapezoidal Triple-GEM detector with wide readout strips was tested in hadron beams at the 
Eermilab Test Beam Eacility in October 2013. The readout strips have a special zigzag geometry and run 
along the radial direction with an azimuthal pitch of 1.37 mrad to measure the azimuthal ^-coordinate of 
incident particles. The zigzag geometry of the readout reduces the required number of electronic channels by 
a factor of three compared to conventional straight readout strips while preserving good angular resolution. 
The average crosstalk between zigzag strips is measured to be an acceptable 5.5%. The detection efficiency 
of the detector is (98.4 ± 0.2)%. When the non-linearity of the zigzag-strip response is corrected with track 
information, the angular resolution is measured to be (193 ± 3) /irad, which corresponds to 14% of the 
angular strip pitch. Multiple Coulomb scattering effects are fully taken into account in the data analysis 
with the help of a stand-alone Geant4 simulation that estimates interpolated track errors. 
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1. Introduction 

Spatial resolution is always an important parameter for position-sensitive particle detectors, such as 
Micro-Pattern Gaseous Detectors (MPGD). It is well known that a Gas Electron Multiplier (GEM) [1] 
detector with straight readout strips of 400 pm pitch can reach 50-70 pm resolution ElE]- In applications 
with large-area GEM detectors, the use of readout strips with such small pitch will quickly lead to a large 
number of required electronic channels and consequently to significant system cost incurred by the readout 
electronics. Readout structures employing pads or strips with a chevron or zigzag geometry were proposed 
and studied in the past to address this issue for MPGDs and other gaseous detectors I1I3E1I71IH]. A zigzag 
strip covers a larger area than a standard straight strip, so that zigzag strips can reduce the number of 
electronic channels needed to read out a given detector area. The challenge then is to make sure that the 
spatial resolution is still adequate and that the crosstalk between strips is well controlled. 

^This draft has been submitted to NIMA for review. 
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We have previously demonstrated for small Triple-GEM detectors that parallel zigzag strips achieve 
good spatial resolution as the sensitivity of the charge sharing among strips to the hit position is enhanced 
between the interleaved “zigs” and “zags” of adjacent strips [9]. Here we expand this approach to large-area 
GEM detectors with a radial readout strip geometry that is appropriate for detector systems with a disk 
or ring geometry as is commonly used in collider experiments. We report results from a beam test of a 
1 -meter-long GEM detector read out with wide radial zigzag strips. 

This paper is organized as follows: In section the design of the zigzag readout board is described in 
detail; in section the beam test configuration and data acquisition system are introduced; in section 
the data analysis methods and basic performances of the zigzag GEM detector are discussed; in section 
the angular resolution, i.e. the spatial resolution in azimuthal ^-coordinate, of the zigzag GEM detector is 
studied under different operating conditions. 

2. The GEM detector and the zigzag readout strips 

The 1-meter-long Triple-GEM detector has the shape of a trapezoid with 22 cm width at the narrow end 
and 45 cm width at the wide end (Eig. [^. It is based on a prototype for the GMS muon GEM upgrade at 
LHG [IQIITT], but has a modified readout board. The gas gaps (drift, transfer 1, transfer 2 , induction) in 
this GEM detector are 3/1/2/1 mm. 

On the printed circuit board (pcb) used for the readout, the zigzag strips are radially divided into eight 
sectors with lengths 10 cm (15 cm) at the narrow (wide) end of the trapezoid and with 128 strips per 
sector (Eig. [^. In a collider geometry these sectors would correspond to different ranges of pseudorapidity 
77 . Erom the narrow end to the wide end, the sectors are labeled from one to eight. The reason for this 
partitioning is to provide some coarse information on the radial position of hits and to have readout strips 
with reasonably low capacitance. In Eig. the details of the zigzag structure are shown. The zigzag strips 
run in radial direction and measure the azimuthal (p coordinate. The opening angle between the first strip 
and the last strip in a sector is 10° and the angular pitch between two neighboring strips is 1.37 mrad. The 
spacing between zigzag tips on adjacent strips is 0.1 mm and the distance between two tips in one strip is 
0.5 mm. These geometric parameters are the same for all strips in all eight sectors on the readout board; 
consequently the strip width increases from ~2.5 mm at the narrow end to ^4.5 mm at the wide end. The 
total number of zigzag strips on the readout board is 1,024, to be compared with, e.g., 3,072 straight strips 
in the standard GMS GEM detector that the zigzag detector is derived from m- This means that 2/3 of 
the electronics can be eliminated by using zigzag strips. 

A miniaturized ceramic high voltage divider is used for powering the GEM detector. The HV divider has a 
resistor chain with a total resistance of 4.4 MO and connects to each electrode of the GEM detector through 
eight pads on the drift board. One HV input channel is connected to the drift electrode and potentials 
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between any two neighboring detector electrodes are provided through resistors in the HV divider. More 
details of the HV divider can be found in reference m- 

3. Beam test configuration and data acquisition system 

The zigzag GEM detector was tested as a tracking detector at the Fermilab Test Beam Facility (FTBF) in 
October 2013. The eRD6-FLYSUB consortium [^conducted this beam test with a variety of GEM detectors. 
Fig. (top) shows a diagram of the detector configurations in the tracking system, where a set of four GEM 
detectors with 2-D readout acts as a reference tracker (noted as “Ref” in the figure). The zigzag GEM 
detector was installed on a movable table in the center between the tracker detectors (Fig. bottom). All 
detectors were operated with an Ar/G02 70:30 gas mixture during the entire test period. 

The RD51 Scalable Readout System (SRS) [12] was used to read out all ten detectors in the beam 
test. With 64 frontend hybrids carrying APV25 chips [13] and four SRS Frontend-Goncentrator/ADG 
combinations, 8,192 channels were read out simultaneously (Fig. [^. Strip signals were digitized at 40 
MHz. Readout was triggered by coincidence signals from plastic scintillators placed in the beam line. Data 
were transferred via Gigabit ethernet to a PC that acquired them using the DATE software and monitored 
them online in the Automatic MOnitoRing Environment (AMORE) (both software suites were originally 
developed for the ALICE experiment and later adapted by the RD51 collaboration for the SRS) [T4] . 

4. Data analysis and offline results 

The data are analyzed with the AMORE software package as well; raw events are decoded and basic 
information on the events are retrieved. Events with multiple hits in any detector are excluded in the 
analysis for simplicity. Beam profiles are checked first to characterize the impinging beams. Then basic 
performance characteristics of the zigzag GEM detector are analyzed, such as strip multiplicity of strip 
clusters, i.e. the number of strips with induced charge on the strips above a certain threshold, strip-cluster 
charge distribution and detection efficiency. For simplicity, we will refer to a strip cluster as a “cluster” 
throughout the remainder of this paper. Crosstalk between zigzag strips is also investigated in this section. 

4 . 1 . Beam profiles 

During the beam test, we tested detectors with the primary 120 GeV/c proton beam and with secondary 
beams containing mixed hadrons (mainly pions with an admixture of kaons and protons) at 20, 25, and 
32 GeV/c momenta. Most of the data for the zigzag GEM detector were taken with mixed hadron beams. 

^eRD6-FLYSUB is a consortium of researchers from FLorida Institute of Technology, Yale University, Stony Brook Uni¬ 
versity, University of Virginia, Brookhaven National Lab, formed to carry out R&D of tracking and particle identification 
detectors for a future electron-ion collider. 
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The 2D hit maps of the tracker detectors (Fig. show that the beam spot of the pure proton beam is 
^2 cm in diameter, while in the case of mixed hadron beams it is ^7 cm in X direction (horizontal) and 
^3 cm in Y direction (vertical). 

4 . 2 . Basic performance characteristics of the zigzag GEM detector 

4 . 2 . 1 . Cluster charge vs. HV 

Before taking data for any scenario, we took a pedestal run with 5k events at standby detector voltages 
so that there was no gas gain. To find a cluster in data runs with full voltage, we set a 5cr threshold cut where 
a is the width of the pedestal distribution. If N contiguous strips are fired in a detector with individual 
strip charges higher than the 5cr cut, then we will consider this as a cluster with strip multiplicity N in the 
detector and accept this event. 

The zizgzag GEM detector was tested at different high voltages Vdrift applied to the drift electrode when 
a 25 GeV/c mixed hadron beam was impinging on central readout sector number five (Fig. 0 . Fig. (top) 
shows the total cluster charge distribution at Vdrift = 3200 V in terms of ADC counts (1 ADC ^ 0.0371 fC). 
The measured cluster charge distribution fits well to a Landau function. Fig. (bottom) shows the most 
probable values (MPV) of the Landau fits as a function of Vdrift; this yields an approximately exponential 
curve as expected for an MPGD. The gas gain during the beam test is estimated to range from 660 (at 
Vdrift = 3000 V) to 10^ (at Vdrift = 3400 V) based on the measured mean cluster charge. After the beam 
test, we measured the gas gain of this detector operated with Ar/C02 70:30 in our lab in Florida (sea level) 
with X rays (AMPTEK Mini-X X-ray generator with gold anode, set at 10 kV and 5 juA) impinging on 
central sector five. A good rate plateau is obtained and the gain reaches close to 3 x 10^ at Vdrift = 3600 V 
(Fig.0 under those conditions. 

4 . 2 . 2 . Strip multiplicity vs. HV 

The strip multiplicity distribution of clusters measured at Vdrift = 3200 V in central sector five is shown 
in Fig. 1^ (top); the mean strip multiplicity is 1.63. This mean value is found to increase quadratically with 
Vdrift (Fig. i bottom), and at the highest tested voltage (3400 V) its maximum is just below three strips. 
This is expected since the strip width in that radial region is ^3.6 mm. 

4 . 2 . 3 . Detection efficiency vs. HV 

The detection efficiency of the zigzag GEM detector for mixed hadrons, which are close to minimally 
ionizing in this test beam, is obtained from the ratio of the number of observed hits to the total number 
of triggered events. The measured efficiency is (98.4 ± 0.2)% on the plateau when a 5cr threshold cut is 
applied to the measured strip charge. Fig. shows efficiencies with different threshold cuts demonstrating 
that the plateau efficiency is not affected significantly by the applied threshold. 
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^.2.^. Response uniformity from position sean data 

The zigzag GEM detector was tested (at Vdrift = 3200 V) with 20 GeV/c mixed hadron beams impinging 
on two different spots in each readout sector, approximately 60 mm apart in the vertical direction, so that 
the response uniformity of the chamber could be examined. Fig. pT| shows the MPVs of Landau fits to cluster 
charge distributions measured in sectors one to seven at these different positions. The sector eight was not 
tested because it could not be reached by moving the stage that the detector was placed on. The response 
varies by ^25% and is lower in sectors six and seven, which are on the wide side of the trapezoidal chamber. 
This non-uniformity is likely caused by a known slight bending of the drift board that was occurred when 
the GEM foils were stretched during assembly. 

4 . 2 . 5 . Crosstalk among zigzag strips 

We were able to investigate the crosstalk between zigzag strips with the beam test data due to a minor 
mistake that was made in the production procedure of the readout pcb, which accidentally connected strips 
number 63 and 127 in each sector to the ground plane. The large input capacitance presented by these 
grounded strips to the amplifier in the APV chip causes large noise as can be seen, for example, in the 
distribution of the pedestal widths for one of the sectors (Fig. [IT]). These plots also show that the adjacent 
strips 62, 64 and 126 are victims of crosstalk as their pedestal widths are slightly higher than the average 
across all other strips. 

We estimate the crosstalk as y^rms^-^.^-^ — rms^^g/rmSaggressor, where rmsavg is the average pedestal 
width observed for strips not affected by crosstalk, rmSaggressor is the pedestal width on the noise source 
strips (63 or 127), and rrnsvictim is the pedestal width on a strip adjacent to the noise source strip. We find 
an average crosstalk among the zigzag strips of (5.5 ± 0.2)% with an rms width of 1.3% (Fig. [T^ . 

5. Angular resolution studies 

5.1. Multiple Coulomb seattering in deteetor materials during beam test 

In the beam test, the total radiation length of the detector materials in the setup composed of ten 
GEM detectors was about 14% (Table [^, which impacts the spatial resolution analysis due to multiple 
Goulomb scattering (MGS) of the tracks in the material. For example, the rms width of the scattering angle 
distribution due to MGS in that amount of materal is estimated to be 147 /irad for 32 GeV/c hadrons using 
the standard MGS formula m- 

A stand-alone Geant4 m simulation was created to study the impact of MGS and to extract the intrinsic 
detector resolution precisely. The beam test setup as simulated in Geant4 is shown in Fig.[^ The material 
distribution from Table is fully implemented in the simulation, while details such as GEM holes and 
readout strip geometries are not. For the Geant4 physics list, FTFP_BERT is used, which includes the 
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Table 1: Material estimate for the tracking detectors in the FNAL beam test. 


Detector 

Gas gaps 

Window mat./ 

Readout mat./ 

Rad. Len. 


[mm] 

thickness [mm] 

thickness [mm] 

[%^o] 

Tracker 1 

3/2/2/2 

Mylar/^ 0.1 

G10 /kapt on / honeycomb 

0.32 

Tracker 2 

3/2/2/2 

Mylar/^ 0.1 

G10 /kapt on / honeycomb 

0.32 

SBS 1 

3/2/2/2 

Al+kapton 

G10 /kapt on / honeycomb 

0.345 

UVA-lm-GEM 

3/2/2/2 

Mylar/^ 0.1 

GlO/kapton/Rohacell foam 

0.42 

FIT-lm-zigzag-GEM 

3/1/2/1 

PCB/3.175 

GlO/3.175 

3.88 

FIT-30cm 

3/2/2/2 

PCB/3.175 

GlO/2.362 

3.42 

FIT-lOcm-1 

3/2/2/2 

Mylar/^ 0.1 

GlO/2.362 

1.5 

FIT-lOcm-2 

3/2/2/2 

Hony comb /3.175 

GlO/2.362 

1.48 

Tracker 3 

3/2/2/2 

Al+kapton 

G10 /kapt on / honeycomb 

0.345 

Tracker 4 

3/2/2/2 

Mylar/^ 0.1 

G10 /kapt on / honeycomb 

0.32 

Ar/COa 

88 mm 



- 0.66 

Air 

^ 3 m 



^ 1 

Total Material 




14% Xo 


MCS model based on Lewis theory. Perpendicular point-like beams start 20 mm in front of the first tracker 
detector (REFl) at position (x,y) = (0 mm, 0 mm) and different beam momenta and particles are simulated. 

We estimate the magnitude of the MCS effect in each detector by running the simulation with perfect 
intrinsic detector resolutions, but with MCS turned on. The simulated hit positions are histogrammed for 
each reference tracker detector and for the zigzag GEM detector. In the simulation, the zigzag GEM can be 
treated as a 2D detector in a Cartesian system. In Eig.[^ the top plot shows the hit position distributions 
in the horizontal X direction for 25 GeV/c pious, while the bottom plot is for 120 GeV/c protons. Beams are 
widened significantly at lower momenta due to MGS and the hit distribution widens as the beam particles 
move downstream. The scattering angles between detectors can be calculated (Fig, The scattering is 
small at the beginning of the beam and it becomes large at the end of the 3-m long tracking system; the 
overall mean scattering angle between tracker detectors one and four is about 97 /xrad with a 66 /irad rms. 
The exclusive residual width for the zigzag GEM detector due to the MGS effect alone is found to be about 


80 jam in Gartesian coordinates (Fig. 16) when simulating 25 GeV/c pions with all tracker detectors set to 
perfect intrinsic resolution. Here ‘exclusive’ means that the hit of the zigzag detector is excluded from the 
track fit. 
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5.2. Alignment of reference trackers 


The first step in the resolution measurement is an alignment of the four small tracking detectors that 
have a Cartesian X-Y strip readout. The trackers are first aligned to each other in Cartesian coordinates. 
The data sample with the highest statistics is used for this alignment, which was taken with 32 GeV/c mixed 
hadron beams. Only events with a single cluster observed in each tracker are used for the alignment. Any 
strip multiplicity is allowed in the clusters including single-strip clusters. 

The first alignment step is to shift each of the four tracking detectors iteratively in the XY-plane to make 
their origins match each other in that plane. The initial shift parameters are mean values from position 
distributions in X and Y coordinates. In each iteration, straight lines are fitted to the hits in X and Y. 
Residuals are histogrammed for each detector and the residual distributions are fitted with a double-Gaussian 
function. Ten percent of the residual mean value of each detector is taken as the shift parameter in the next 
iteration to avoid over correct ions. The resulting residual mean values converge quickly towards zero after 
40 iterations, as can be seen in Fig. This provides a first coarse alignment. In a second alignment step, 
we correct also for relative rotations of the tracking detectors around the beam in the XY-plane. We again 
fit straight lines to the hits in X and Y and iterate through a succession of offsets and rotations around the 
beam axis relative to the first tracking detector until the residual means from the track fits are very close 
to zero and the of Iho track fits are minimized. In each iteration, the detectors are first shifted and then 
rotated; then new residuals and rotation angles are calculated. After about 20 iterations the residual mean 
values become flat within ±0.4 /am around zero (Fig. [T^. From this step, we get eight shift parameters 
and three relative rotation angles as refined alignment parameters for the tracker system. The final step is 
to optimize the three rotation angles one-by-one. With the shift parameters and two rotation angles kept 
fixed, the third rotation angle is changed in a small range around the value from the second step with a 
1 mrad step. We find parabolic curves of mean of tracks vs. angle and calculate the optimized rotation 


angle from the minimum of parabolic fits to those curves (Fig. 19). The final rotation angle is taken as the 
average of these angles obtained for X and Y directions. 


5.3. Spatial resolution of reference trackers 

The Geant4 simulation of pious traversing the setup allows us to estimate track errors and to measure 
intrinsic detector resolutions when MGS effects are included. The resolutions of the tracker detectors are 
studied to ensure that their resolutions are good enough to serve as precision reference detectors. In the 
first simulation step, the intrinsic detector resolutions are first set to zero (perfect resolution) but with 
MGS turned on, and then the obtained hit positions are additionally smeared by hand using a Gaussian 
to simulate the intrinsic detector resolution. We initially assume that the four tracker detectors have the 
same resolution and smear the simulated hit data with resolutions from 50-80 /am in 5 /am steps. For each 
smeared resolution, we calculate the exclusive track-hit residuals for each detector and get the residual 
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widths, so we can compare exclusive residual widths from simulation with those observed in experimental 
data (Fig. [^. When the residual width matches for a tracker detector, the corresponding input resolution 
used for the Gaussian smearing is taken as the intrinsic resolution of that detector. We do this in both X 
and Y coordinates and the average of the two is taken as the final resolution for a tracker detector. The 
resolutions of the tracking detectors 1-4 are found to be 73, 70, 59, and 68 /im, respectively. The statistical 
uncertainties of these resolutions are smaller than 2.5 fim (half of the step size). These measured resolutions 
fall into a typical range for standard GEM detectors with 400 jam strip pitch. 


5.4- Alignment of GEM detector with zigzag strips 

Since the zigzag strips in the large GEM detector run in radial direction and measure only the azimuthal 
^-coordinate, we need to perform the resolution study for this detector in polar coordinates. The tracker 
hit positions are transferred to the polar coordinate system that is naturally given by the geometry of the 
large zigzag detector with the vertex of the trapezoid used as the origin of this polar system (Eig. [21] ), and 
the tracker tracks are then refit in these polar coordinates. 

In the polar coordinate system, the ^-coordinate of the center of each zigzag strip in the large trapezoidal 
detector is simply given by = —0.5 • a n • (a/(128 — 1), where a = 10*^ is the opening angle of the 
trapezoidal shape and n = 0,1,.. .,127 is the strip number. The ^-position of a hit is initially determined 
from the barycenter, or centroid, of the cluster using the strip charges as weights. 

Eor transforming the Gartesian tracker coordinates into the polar coordinate system, we need to find the 


correct Xoffset and Yoffset from the new origin to the center of the tracker (Eig. 21) using the tracks in the 
^-coordinate. This has to be done for each test scenario. Note that since the tracker is already internally 
aligned, the centers of the individual trackers already match, so there is only a single pair (Xoffset, Yoffset) 
of offsets to be found for each scenario. The idea behind our alignment procedure is that any misalignment 
will shift the residual means away from zero and will increase the residual widths when fitting tracks in the 
^-coordinate with hits in the trapezoidal GEM. 

We use a two-fold iteration loop to find the X and Y offsets for each data-taking scenario. Eirst, we keep 
a fixed Xoffset value and step through Yoffset values with a fixed step of 0.2 mm within a reasonable physical 
range, then change to another Xoffset value with a step of 1 mm and vary the Yoffset again until all Xoffset 
values have been covered. Eor each (Xoffset, Yoffset) pair, tracks are fitted in the ^-coordinate including and 
excluding the 0-hit in the zigzag detector and the corresponding inclusive and exclusive residuals for the 
zigzag GEM detector are histogrammed. We record residual means and widths, as well as the mean of 
the track fits in the ^-coordinate. 

We use a three-step procedure to find the best estimate of the offsets. Eig. [22] shows an example of this 
procedure for the data recorded in ? 7 -sector five of the zigzag GEM at 3300 V. Eor a given Xoffset value we 
record that Yoffset value which produces a well-centered ^-residual for the zigzag detector with a mean of 
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zero (Fig. left). Then, for a given Yoffset? we plot the residual widths as function of Xoffset and find the 
Xoffset value that minimizes the residual width using a parabolic fit (Fig.center). This yields two sets of 
(Xoffset: Yoffset) paus that are plotted as two curves. The (Xoffset: ^offset) values for which those two curves 
intersect is taken as the best estimate of the alignment offsets (Fig. right). Best (Xoffset: Yoffset) pairs 
based on track-y^ are also plotted and confirm the results. 

Finally, it is not guaranteed that tracker and zigzag detector are installed without a relative rotation in 


the X-Y plane (Fig. 23 , top). Since aligning the zigzag GEM detector also requires matching the X direction 
to the trackers so that all calculations of (j) are correct, the rotation angle of the zigzag GEM detector needs 
to be checked. The mean of tracks in (j) vs. rotation angles is a parabola and the angle can be calculated 
from the minimum; it turns out to be very close to zero (Eig. bottom). 


5 . 5 . Corrections for non-linear strip response of the GEM detector with zigzag strips 

Before proceeding to the measurement of the angular resolution for the zigzag GEM detector, we correct 
its hit data. The reason to do so is that there is a non-linear response of the zigzag strips when the hit 
positions are calculated from the centroids of clusters. Specifically, we use track information to apply the 
non-linear corrections to the hit positions. Eor each cluster with a strip multiplicity N > 1 , we define a 
quantity 77 = Sc — Smax: where Sc = • Si/E^^^^qi is the centroid position of the cluster in terms of 

strip number; si and qi are the number and the induced charge (in ADG counts) for the strip in the 
cluster and Smax is the number of the strip in the cluster on which the maximum charge is induced. The 
quantity 77 is then a measure of the difference between centroid position and the maximum-charge strip 
independent of the absolute strip number. We treat different strip multiplicities separately because 77 has 
different characteristics for odd and even strip multiplicities. Eor even multiplicities, Sc tends to be between 
two strips, whereas for odd multiplicities it is close to the center of the strip with maximum charge. In 
practice, we only need to correct clusters with strip multiplicities N = 2 and N = 3 using corresponding 
quantities 772 and 773. Gorrections for N > 3 are not necessary since there are very few of such clusters (see 
also Eig. [^. 

Using reference tracks given by the tracker, exclusive residuals of hits in the zigzag GEM are then plotted 
vs. 77. We observe that the resulting distribution is not flat, which indicates a non-linear response of the 
zigzag strips to the true hit position. As shown in Eig. exclusive residuals are plotted vs. 772 and 773 and 
their profiles are fitted separately to appropriate fit functions. A 10^^-degree polynomial produces a good 
fit for the 772 distribution (Eig. top), while a serpentine function is used for the 773 distribution (Eig. 
bottom). We then globally correct the orginal hit positions by subtracting the 77-dependent offsets obtained 
from the fit functions from the original hit positions. Eig. shows that the exclusive residual vs. 77 plot 
is much flatter after this correction. Gonsequently, the overall residual distribution becomes significantly 
narrower and the spatial resolution is improved when the corrected hit positions are used. Note that the 
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tracker detectors do not need to be corrected since their resolutions are already good enough as they have 
straight strips with much smaller pitch and consequently much more linear responses. For the position scan 
data, we apply individual correction functions for each sector since the strip width changes along the radius. 
For the HV scan data, which were taken in a single sector, the correction functions are obtained from all 
HV scan data combined together and applied equally to different voltage points. We have checked that the 
difference is very small if we find correction functions using only data taken at one voltage point. 


5.6. Angular resolution measurement for GEM deteetor with zigzag strips 

In polar coordinates, we smear the simulated tracker hits with the realistic intrinsic tracker resolutions 
obtained above and with MCS fully taken into account. We fit a linear track to the azimuthal coordinates of 
the four smeared hits from the tracker detectors and simulate an unsmeared hit for the zigzag GEM detector. 
The resulting width of the exclusive residual ^interpoi.track ~ ^unsmear.zigzag distribution for the zigzag GEM 
is then a measure of the interpolated-track error (IE) at the position of the zigzag GEM with MGS taken 
into account. Eor the experimental data, we also fit linear tracks to the ^-coordinate and calculate the 
exclusive residual width (ER) for the zigzag GEM detector. The intrinsic resolution a for the zigzag GEM 
can then be calculated by subtracting the two quantities in quadrature: a = \/ER^ — lE^. We check that 
this method gives accurate results using the simulation by also smearing the resolution of the probed zigzag 
GEM detector from 10 to 390 /irad in 10 /irad steps and by then applying the same method. We find almost 
perfect agreement between input resolution from Gaussian smearing and calculated intrinsic resolution a 
over the full range (Eig. [2^. 

Eig.[^ shows the resulting measured intrinsic angular resolutions for the zigzag GEM detector in central 
sector five at different Vdrift values. The resolution is around 180 /irad on the efficiency plateau for 2-strip 
and 3-strip clusters when the non-linear responses of zigzag strips are corrected. In the bottom right plot, 
we see that the overall angular resolution for the zigzag GEM including single-strip clusters, which cannot 
be corrected, is a = (193 ± 3) /irad at highest tested voltage. This value corresponds to 14% of the angular 
pitch of the radial strips. The resolutions are also measured for different positions on the zigzag GEM 


operated at 3200 V (Eig. 28). 


6. Summary and Conclusion 

The large-area zigzag GEM detector performed very well in the Eermilab beam test. The ^5.5% average 
crosstalk between the zigzag strips is small and does not show significant impact on signal collection. The 
detection efficiency is above 98% for charged particles as expected for a GEM detector. The zigzag GEM 
detector achieves an angular spatial resolution of 193 /irad or ^14% of the angular strip pitch after the 
centroid positions of clusters are corrected for non-linear strip response. 
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In conclusion, the zigzag readout design is considerably more cost-effective for large-area GEM detectors 
than conventional straight strip readout structures while preserving good performance. Consequently, GEM 
detectors with zigzag readout strips are a viable option for application as affordable tracking detectors at 
future colliders such as the Electron-Ion Collider (EIC). 
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Figure 1: The 1-meter-long trapezoidal GEM detector with zigzag readout strips. The detector has eight sectors; each comprises 
128 radial zigzag strips and is read out with an APV25 hybrid board. 



Figure 2: A photo of the side of the readout pcb that features the radial zigzag strips. The inset is an image of the zigzag 
structure taken with a microscope. The strips run in radial direction with an angular pitch of 1.37 mrad. The distance between 
two tips in neighboring strips is 0.1 mm, and the distance between two tips in the same strip is 0.5 mm. 
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Figure 3: Top: The overall configuration of all tracking GEM detectors in the beam line (not to scale). Bottom: A photo of 
the zigzag GEM detector (pointed out by the yellow arrow) in the beam line. 
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Figure 4: The Scalable Readout System used in the beam test. Top: HDMI cables from APV25 hybrids connect to ADCs that 
are connected to Front-End Concentrators (FECs). Bottom: The FECs in turn communicate with a DAQ PC through Gigabit 
ethernet. One ADC+FEC unit reads out 16 APV25 hybrids, which corresponds to 2,048 channels. 
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(c) 32GeV/c mixed beam (d) 120GeV/c proton beam 

Figure 5: 2D maps of hits recorded in the first tracker detector for different beam momenta. Only events with hits in all 
four trackers are shown. As expected, the beam spot of the primary proton beam is much smaller than the beam spot of the 
secondary mixed hadron beams. The horizontal and vertical lines visible on these plots are due to grid spacers in the third 
tracker detector that locally reduce the detection efficiency. 
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Figure 6: Top: Cluster charge distribution measured in central sector 5 of the zigzag GEM detector at V^rift = 3200 V with 
25 GeV/c mixed hadron beam and fitted to a Landau distribution. Bottom: Most probable values (MPV) of the measured 
Landau distributions vs. V^rift- (1 ADC count corresponds to about 0.0371 fC.) 
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Figure 7: Rate and gain vs. V^rift measured in central sector 5 of the zigzag GEM detector with Ar/C02 70:30. As opposed 
to all other measurements discussed in this paper, this measurement was performed after the beam test at sea level in our lab 
in Florida with x-rays. The x-rays were generated by an AMPTEK Mini-X x-ray gun with gold anode that was operated at 
10 kV accelerating voltage. Errors are smaller than marker sizes. 
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Figure 8: Top: Strip multiplicity for clusters in the zigzag GEM detector measured at Vdrift = 3200 V in central sector 5 with 
a 25 GeV/c mixed hadron beam. Bottom: The mean strip multiplicity vs. Vdrift follows a quadratic function. Errors are 
smaller than marker sizes. 
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Figure 9: Detection efficiency in central sector 5 of the zigzag GEM detector for 25 GeV/c mixed hadron beam as a function of 
Vdrift s-iid for different hit thresholds (errors are smaller than marker size). The applied thresholds are multiples of the width 
(cr) of the pedestal distributions for individual strips. The data are fit to sigmoid functions as given in the inset to determine 
the efficiencies on plateau. 
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Figure 10: Most probable value of cluster charge distributions from Landau fits in different sectors (Vdrift = 3200 V). Two 
points were measured in each sector. “Middle” means the position is roughly in the center of the sector and ~60 mm lower 
than the “Upper” position in that sector. Errors are smaller than marker sizes for higher sector numbers. 
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Figure 11: Pedestal widths (rms) observed for zigzag strips in the vicinity of two noisy strips (“aggressors”) in one sector. 
Pedestal widths in adjacent strips (“crosstalk victims”) are slightly higher than the average pedestal width due to crosstalk 
from the aggressors. 
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Figure 12: Overall crosstalk measured for crosstalk victims in two pedestal runs and in several sectors of zigzag GEM. 



Figure 13: The stand-alone Geant4 geometry of all GEM detectors in the FNAL 2013 beam test. The large-area GEM detector 
with radial zigzag strips is labelled as FITGEM. 
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Figure 14: Simulated hit positions in the tracker and in the zigzag GEM detector (FITGEM) for 25 GeV/c pions (top) and 
120 GeV/c protons (bottom). 
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Figure 15: Distributions of scattering angles between different reference detectors due to multiple scattering. The mean and 
rms values in firad of each distribution are indicated on the plot. 



Figure 16: Simulated exclusive residual distribution for the zigzag GEM detector for 25 GeV/c pions. All detectors are set to 
have perfect intrinsic resolution here so that the width of this distribution is a measure of the MGS effect at the position of 
the zigzag GEM. 
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Figure 17: Residual means vs. iteration number in the first step of aligning the tracker. In this step, the detectors are only 
shifted in each iteration. The inset is a zoom-in around zero. 



Figure 18: Residual means vs. iteration number in the second step of aligning the tracker. In this step, detectors are alternately 
shifted and rotated relative to tracker detector REFl in each iteration. 
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Figure 19: Track-x^ in X-Z (top) and Y-Z (bottom) planes vs. rotation angle of tracker detector REF2 relative to tracker 
detector REFl. This is an example for the third step of aligning the tracker; the final relative rotation angle obtained from 
averaging the minima of the two parabolic curves is 5.72 mrad. 
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Figure 20: Exclusive tracker detector residuals for X (left) and Y (center) coordinates obtained with common input values for 
the intrinsic detector resolutions (Gaussian smearing in 5 jim steps, black points) in the Geant4 simulation with MGS compared 
with experimental residuals (blue points) in the beam test. The numbers in blue (left and center) are resolutions for the tracker 
detectors in the beam tests obtained from this comparison with the simulation. When feeding these resolutions (averaged over 
X and Y) for each tracker detector back into the simulation as inputs for smearing, the resulting simulated residuals are found 
to be consistent with the experimental data (right) 



Figure 21: Schematic diagram for transformation of the Gartesian tracker coordinates into the natural polar coordinate system 
of the zigzag GEM detector, which has the vertex of the trapezoid in its origin. 
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Figure 22: Optimization of X and Y offsets for tracker detectors for beam impinging on zigzag GEM in central sector 5. Left: 
At fixed Xoffset = —1876 mm, the residual mean vs. Yoffset curve is a line; the best Yoffset cit this Xoffset point is calculated 
from requiring that the residual mean be equal to zero: Yoffset = —35.05 mm. Center: At fixed Yoffset = “35 mm, the residual 


width vs. Xoffset is a parabola; the best Xoffset is calculated from the minimum of this parabola: Xoffset = —1876.6 mm. Right: 
Scatter plots of best (Xoffset 5 ^offset) points from residual mean vs. Yoffset and residual sigma vs. Xoffset 5 as well as plots of 
vs. Xoffset and vs. Yoffset 5 which strongly overlap with the former. 




Figure 23: Top: Schematic diagram for slight relative rotation of the zigzag GEM detector. Bottom: Mean x^ of tracks in cj) 
coordinate vs. rotation angles for the same data as in Fig. |22| from the minimum of the parabola the global rotation angle of 
the zigzag GEM is found to be —51 /xrad. 
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Figure 24: Top: Exclusive residuals of 2-strip clusters in zigzag GEM vs. 772 fitted with 10-degree polynomials. Bottom: 
Exclusive residuals of 3-strip clusters in zigzag GEM vs. 773 fitted with a serpentine function. 


27 


































3 0 
^• 0.0002 
-0.0004 
-0.0006 r 

- 0.0003 ‘ ^ 

ft ftft4 Fi I L I I ] I L I I iV I I I I I ] I [ I ab I I I I I I I I t I 1 > L I I I I 1*1 I I I j I L i‘J ^Mn 

-6,5 -0.4 -0.3 -0.2 -0.1 0 0,1 0.2 0,3 0,4 0.5 

n of 2 and 3-Strip clusters 


0.001 


feslduaJ ela attef | 

Enirj«£ 

14925 1 

Maarix 

(1.D01642 

Maany 

7.97Se-06 


02544 





- 0.0002 

■0.0004 


■0.0006 
-0.0008 r 


■ ■ ■ . ■ ■ ■ ' I 


E 

10 u 

z 


ft ftft-j Lli I J t I J 1 L I I I 1 I I I I J I I! I I L I I i L I j‘l L I I I t I I I I J I Li J J L L J a L I 

-6,5 -0,4 -0.3 -0,2 -0,1 0 0.1 0.2 0,3 0,4 0.5 

11 Of 2 and 3-strip clusters 

Figure 25: Exclusive residuals of zigzag GEM vs. r] for strip multiplicities N=2 and N=3 clusters before (top) and after (bottom) 
correction. 



Figure 26: Intrinsic angular resolutions calculated using track error estimate for the large GEM detectors vs. the smeared input 
resolutions in simulation, showing very good agreement. Errors are smaller than marker size for most of the points. 
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Figure 27: Angular resolutions vs. V^rift measured in central sector 5 of the zigzag GEM before (black) and after (blue) 
non-linear response corrections. The plots show results for different strip multiplicities of the clusters. 



HV M 


ouu 

280 
260 
240 
220 
200 
180 
160 
H An 


r- 

With 2 and 3-strip clusters 

▼ Before correction 

♦ After correction 













- 







T 

































.j 

f. 





f 

. i 

i . 

I 





. j 

\ . 



f . 

:_I 

i_i 

i 

i_i 

f . 

I_I 

i_I 

i_I 

i_i 

i_i 

I_I 

i_ 


29 


































































































































































Figure 28: Angular resolutions measured at V^rift = 3200 V in different positions of the zigzag GEM before and after corrections 
using all strip multiplicities. The top (bottom) plot is for positions in the central (upper, i.e. ~60 mm above central) part of 
each sector. 
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